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Summary 

Finite-element, transient heat transfer analyses were 
performed for the first-stage blades of the space shuttle main 
engine (SSME) high-pressure fuel turbopump. The analyses 
were based on test engine data provided by Rocketdyne. Heat 
transfer coefficients were predicted by performing a boundary- 
layer analysis at steady-state conditions with the STAN5 
boundary-layer code. Two different peak- temperature over- 
shoots were evaluated for the startup transient. Cutoff transient 
conditions were also analyzed. A reduced gas temperature 
profile based on actual thermocouple data was also considered. 
Transient heat transfer analyses were conducted with the 
marc finite-element computer code. 

Introduction 

The purposes of this study were to calculate the thermal 
response in an SSME turbine blade for subsequent structural 
analysis and to gain greater understanding of blade damage 
modes, convective cooling effects, and the effects of changes 
in gas profile. The heat transfer analyses were performed by 
using a nonlinear, three-dimensional, finite-element analysis 
computer code. The calculated metal temperature distribution 
and history were then used for a nonlinear structural analysis 
of the blade with the same three-dimensional, finite-element 
analysis computer code. 

Experimental measurements of gas temperature profiles in 
the ssme turbopumps are difficult to obtain because of the 
severe operating environment within the SSME. However, the 
operating temperatures of the turbine blade in rocket engine 
turbopumps are primarily a function of the hot gas flow and 
cooling. The temperature field is determined by the heat 
transfer from the hot gas to the blade. This heat transfer and 
its variations are determined from knowledge of the gas film 
coefficients. The temperature gradient through a blade is 
governed by such factors as blade geometry; location of 
shrouds, platforms, dampers, etc. ; and the presence of coating 
material. Determining blade temperatures at the startup and 
cutoff of an engine cycle is difficult because of the complex 
nature of flow through an accelerating turbine. Thus, the time- 
temperature history profile at the startup transient, during 
steady state, and at cutoff was obtained through a combination 
of analytical and experimental results. 


In the past, first-stage turbine blades in the high-pressure 
fuel turbopump (hpftp) of the space shuttle main engine have 
undergone cracking near the blade shank region and at the 
airfoil leading edge adjacent to the platform. These cracks 
apparently were initiated during the first few mission cycles 
by severe thermal transients and were propagated by vibratory 
excitation. To achieve the necessary durability, these blades 
are currently being cast by using directional solidification. 
Single-crystal alloys are being considered for future SSME 
applications. Since the primary interest herein was in the airfoil 
region, the platform and shank were excluded from the finite- 
element model to reduce the size of the problem and the related 
computing time. Also, knowledge of flow conditions at the 
platform, shank, and fir tree was lacking. The flow conditions 
at the shank and fir tree involve the mixing of hot main-stage 
gas with hydrogen coolant flowing up the disk. Details of this 
process are not known, and hardware evidence indicates that 
the thermal environment in these regions varies widely. 

First-stage blades in modem, high-pressure gas turbines 
generally require film cooling. However, adding a layer of 
cooling air over the surface reduces turbine efficiency. To 
minimize this loss in efficiency, it is necessary to predict the 
external heat transfer coefficient during the airfoil design. 
Stepka (ref. 1) shows that the greatest improvement in the 
accuracy of metal temperature, and hence blade life prediction, 
can be achieved by improving the prediction of the external 
heat transfer coefficients. Therefore, accurate boundary-layer 
analysis and well-defined surface heat loads, aerodynamic 
losses, and boundary conditions will substantially improve the 
accuracy of heat transfer analyses. This report addresses the 
problems of calculating the thermal response of high- 
temperature gas-path components, such as turbopump blades 
in space power propulsion systems. The first-stage HPFTP 
blade in the SSME was selected for this study. A three- 
dimensional, finite-element heat transfer analysis based on an 
engine test cycle was performed by using the MARC finite- 
element code (ref. 2). 


Analytical Procedure 

A three-dimensional, finite-element model of the airfoil 
(fig. 1) was constructed of eight-node, solid, isoparametric 
elements. This model, consisting of 360 elements with 576 


l 



TRAI LING-EDGE BASE —y 



^ LEADING-EDGE BASE 
Figure 1.— Airfoil finite-element model. 

nodes, was a shortened version of a finite-element model 
created by Lockheed for a NASTRAN steady-state structural 
blade analysis (ref. 3). The main difference between the NASA 
Lewis and Lockheed models was that the blade base and most 
of the platform were omitted for the MARC heat transfer 
analysis in order to reduce the size of the problem as well as 
the computing cost for subsequent cyclic structural analysis. 
Analysis of the shank region was also hindered by lack of 
knowledge as to its thermal environment. The airfoil had a 
span length of 2.2 cm and a span-to-chord- width aspect ratio 
of approximately unity. 

The finite-element formulation in MARC for the temperature 
at time t can be expressed as 


T(x) = N(x) T (1) 


Temperature T_within the element is interpolated from the 
nodal values T of the element through the interpolation 
functions N as given in equation (1). The governing equation 
of the heat transfer problem is 


C(T)T + K(T)T=Q (2) 


In equation (2), C( T) and K(T) are the temperature-dependent 
heat capacity and thermal conductivity matrices, respectively; 
T is the nodal temperature vector; t is the time derivative of 
the temperature vector; and Q is the heat flux vector. 

The backward difference scheme selected for discretizing 
the time variable in equation (2) yields the following 
expression: 

C(T)T n -x ( 3 ) 

Equation (3) computes nodal temperatures for each time-step 
increment. 

For evaluating temperature-dependent matrices, the tem- 
peratures at two previous steps provide a linear (extrapolated) 
temperature description over the desired interval. 

T(t) = T(t - At) = — T T(t - At) - T(t - 2Af)l (4) 


This temperature is then used to obtain average thermal 
properties of the material /over the interval to be used in 
equation (3). 


-J-C(D +tf(D T n = Q„ + ± 

At J At 



( 5 ) 


During each iteration, the average properties are obtained from 
the results of the previous iteration. 

A boundary-layer analysis was conducted to compute the 
gas film coefficients needed for the heat transfer computation. 
The first step in the analysis was to use the inviscid flow code 
MERIDL (refs. 4 and 5) to determine the flow characteristics 
at the plane midway between the two mean-camber lines of 
adjacent blades. The second step was to use the tsonic code 
(ref. 6) to calculate each of the 1 1 streamlines blade to blade. 
By using the meridl and TSONIC codes, a quasi-three- 
dimensional calculation was effected. 

The thermal environment of the blade was determined by 
employing a two-equation Prandtl mixing-length turbulence 
model to conduct boundary-layer calculations throughout the 
flow. The modified STAN5 code (ref. 7) was used to compute 
the boundary-layer equations as well as the heat transfer 
coefficients. The theoretical basis of this computation is 
described in detail in reference 8. The blade was assumed to 
operate at the steady-state portion of a full power cycle. The 
mission cycle (fig. 2) consisted of a 4.5-sec transient from 
startup to steady state, a 23.5-sec hold at steady state and a 
4-sec transient cutoff or shutdown. Rocketdyne provided the 
gas conditions from engine test data. 


2 


2000 | — 



r UPSTREAM BOUNDARY DOWNSTREAM BOUNDARY 

\ OF ORTHOGONAL MESH OF ORTHOGONAL MESH ~J 

\ AIRFOIL TIP / 



Z DIRECTION 

Figure 3.— Orthogonal mesh of HPFTP first-stage-blade airfoil. 


The finite-difference flow mesh used in the computation is 
shown in figure 3. The airfoil tip and hub cross-sectional 
geometries that were input into the flow codes are shown in 
figure 4. Rocketdyne supplied the inlet gas temperature and 
pressure distributions (table I). The STAN5 program uses air 
properties. To simulate the hydrogen/oxygen flow phenom- 
enon, the Reynolds number was adjusted to maintain dynamic 
similarity. The leading edge of an airfoil poses a special 
problem for the STAN5 program in that a boundary-layer 
profile is required and therefore the calculations cannot start 
from a stagnation point. Gaugler (ref. 7) modified the STAN5 
program by adding a subroutine to calculate the laminar 
boundary-layer profile on a cylinder at a specified surface 
distance from the stagnation line. However, when an attempt 
was made to get detailed heat transfer results in the cylindrical 
leading-edge region, the heat transfer coefficient predicted in 
that region did not agree with experimental correlation. Up 
to about 30° from stagnation, the STAN5 results were strongly 
dependent on the starting point, 0 O . 

Therefore, to compute heat transfer coefficients for the 
leading edge, Gaugler recommends the use of the Martinelli 
et al. (ref. 9) correlation for local heat transfer on a cylinder 
in a laminar crossflow as given by 


Nu d =1.14 Re d 5 Pr 0A 



for <f> < 80° (6) 



Figure 4.— Cross-sectional geometries of airfoil tip and hub sections. 

heat transfer for a cylinder in crossflow was used. It is 
expressed by 


Nu = 0.3 + 0.62 Re 05 Pr 033 




where Nu d is the Nusselt number based on the leading-edge 
diameter, Re d is the Reynolds number based on the free- 
stream approach velocity and the leading-edge diameter, Pr 
is the fluid Prandtl number, and 0 is the angle, in degrees, 
from the stagnation point. 

Equation (6) is usually applicable only to laminar flow 
calculations for high-speed gas turbines. However, it was 
introduced into this analysis since angle effects were included. 
To account for transitional effects into turbulence, the average 


Equation (7) was recommended by Churchill and Bernstein 
(ref. 10). This correlation appears to provide a lower bound 
for Re Pr > 0.4 and a reasonable approximation for all Re 
and Pr. It differs from the prior correlating equation for forced 
convection to cylinders in that it provides a varying power 
for Re and Pr. The effective power of Re in equation (7) 
depends on Pr. Note that Nu is the average for the cylinder. 

Calculations of the local values for the heat transfer coeffi- 
cients at the leading-edge stagnation region were implemented 
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TABLE I.— ROCKWELL INTERNATIONAL SSME HPFTP TURBINE 
GAS-PATH CONDITIONS 

[Reference, 12-7-83 balance at full power level. Shaft speed, 36 742 rpm; mean 
blade speed, 498 m/sec; oxygen/fuel mixture ratio, 1.0364; gas constant, 625.2 
J/kg K; viscosity, 2.852 x 10~ 5 Pa sec; pressure ratio, 1.4491.] 


Condition 

Stage 

First 

Second 

Blade row 

Nozzle 

Rotor 

Nozzle 

Rotor 

Location 

1 

2 

3 

4 

5 

Pressure, MPa: 







Total 

38 


31.52 


26.29 

Static 


33.47 

30.94 

27.56 

25.76 

Temperature, °C: 







Total 

832 

774 


786 


743 

Static 


798 


780 

752 

738 

Specific volume, m^/kg 


0.068 

0.072 

0.079 

0.0828 

Specific heat at constant 

7.86 

7.82 

7.79 

7.77 

7.75 

pressure, kj/kg K 







Ratio of specific heats 

1.3467 

1.3492 

1.3505 

1.3524 

1.3534 

Compressibility 

1.0587 

1.0526 

1.0491 

1.0444 

1.0419 

Pressure ratio 



1.2086 


1.1990 

Reaction, percent 



37.1 


36.6 

Work, kJ 




386 


345.63 

Power, kW (percent) 


29 961 (52.8) 

26 824 (47.2) 


by multiplying the results obtained from equation (6) by the 
amplification term represented in equation (8). The 
multiplication combined the effects of angle, turbulence, and 
turbulent flow 


A=f 


' Re Y’ 62 T 
282 000/ 


( 8 ) 


where/ 5 accounts for turbulence factor and the rest of equation 
(8) for turbulent flow. Nominal turbulence intensity Tu values 
for rocket engines are 10 percent or higher according to 
reference 11. Because of the high-temperature excursions 
associated with the SSME environment, the Tu for the present 
work was chosen to be 10 percent. Therefore, on the basis 
of an Nu/Re 0 5 equal to 1.1286 Pr 0A and Tu Re 0,5 equal to 
48 (at a Tu of 10 percent), a value of 1.6 was approximated 
for / from figures 7 and 8 of reference 1 1 . 

The calculated heat transfer coefficients along the airfoil 
surfaces from the leading to the trailing edge are shown in 
figure 5. The heat transfer coefficients varied by no more than 
20 percent at any axial location over most of the airfoil surface. 
Boundary-layer analysis showed higher heat transfer 


coefficients at the leading edge than at the trailing edge. 
Because of flow compression at the stagnation point, the heat 
transfer coefficients varied primarily from rapid changes in 
surface velocity. 

To simplify the airfoil heat transfer analysis, the following 
assumptions were made: (1) radial hot-gas temperature profiles 
were very small, (2) circumferential variations in hot-gas 
temperature distribution were also small, and (3) the blade tip 
was insulated. Boundary conditions that reflected cooling by 
hydrogen fuel at the blade-to-disk attachment region were 
imposed by activating the MARC user subroutine FORCDT. The 
gas temperature was assumed to be constant over the airfoil 
surface at each time step. Boundary-layer predictions at full 
power level conditions were input into the thermal model by 
means of the marc user subroutine film. Similarly, hot-gas 
temperature and Reynolds number histories were incorporated 
in the input. Figures 6 and 7 reflect 109 percent rated power 
level conditions. 




SURFACE DISTANCE FROM LEADING EDGE, X, MM 

(a) At airfoil root. 

(b) At airfoil tip. 

Figure 5.— Distribution of heat transfer coefficients. 
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Figure 6.— Turbine inlet gas temperature (cutoff transient). 

To obtain transient boundary conditions for both heated and 
cooled surfaces along the airfoil, these boundary-condition 
values were scaled according to the transient flow and 
temperature. In this scaling procedure, the flow was assumed 
to be turbulent. The steady-state heat transfer coefficients were 
adjusted by a factor that is the ratio of the Reynolds number 
over the steady-state Reynolds number. The Reynolds number 
ratio was raised to the power of 0.8, a value assumed on the 
basis that correlations for heat transfer coefficients along a 
flat plate (which is the basis of this analysis) indicate that the 
heat transfer coefficient is proportional to the power of 0.8 
for turbulent flow. In designing film-cooled gas turbine blades, 
this procedure has been found to be appropriate and yields 
satisfactory results. 

Two different peak- temperature spikes, or overshoots, were 
evaluated for the startup transient. (These peak transients were 
based on data that Rocketdyne supplied for the existing SSME. 
The blades were cast from directionally solidified mar-m 246 
alloy. Temperature-dependent properties for this alloy were 
provided primarily by Rocketdyne.) The first spike was caused 
by the liquid-oxygen-rich mixture that occurred when the 
sudden increase in fuel prebumer chamber pressure at ignition 
reduced the fuel flow. The second spike was caused by fuel 
oscillation, which is a cyclic phenomenon observed in the 
SSME. This oscillation reduced the fuel flow to the fuel 
prebumer at about 1.3 sec. As the engine pressure increased, 
the oscillation vanished and the engine reached the main-stage 
operation. The notable feature of the total temperature profile 
for the maximum temperature spike condition (fig. 2) is the 
temperature excursion to 1800 °C at 0.6 sec. Because the 
thermocouple lagged in its response to a fast-changing transient 
gas temperature, this profile is a result of corrected thermo- 
couple measurements and not the actual data readings. The 
correction brought the gas temperature up to 1800 °C. A 
reduced temperature profile was also evaluated in which the 
same spike was limited to 900 °C. 

In this analysis, the automatic time-stepping scheme in 
MARC was used. It was based on a maximum temperature 
change per increment of 10 deg C. The program adjusts the 


5x10 5 



2* 1 i I 

29 30 31 32 

TIME, sec 


(a) At startup transient. 

(b) At cutoff transient. 

Figure 7.— Reynolds number. 

time steps to conform to this criterion according to the scheme 
defined in volume F of the MARC user manual (ref. 2). 
Tolerances are placed on the maximum temperature change 
before the program recalculates nonlinear effects and on the 
maximum temperature variation between the temperature used 
to evaluate properties and the resulting solution. The automatic 
time stepping resulted in 180 increments, and it took 755 sec 
of central processing unit time to run the problem on the Lewis 
Cray xmp computer. 

Discussion of Analytical Results 

Transient metal temperatures were obtained for the two 
startup spikes previously defined. The transient metal 
temperature responses at the hot-spot locations (the leading- 
edge base, the trailing-edge base, and the trailing-edge tip 
(fig. 1)) are plotted in figure 8. Figure 9 shows the response 
of the same hot-spot locations for the reduced temperature 
spike. 

At 0.6 sec after startup (fig. 10(a)), when the first ignition 
spike evaluated in the analysis occurred, the temperature 
gradient was essentially chordwise. This gradient tended to 
decrease until it reached the leading edge, where a significant 
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(a) At startup (corrected gas temperature profile), 
(b) At cutoff. 

Figure 8.— Airfoil transient temperature. 



Figure 9.— Airfoil transient temperature at startup (uncorrected gas temper- 
ature profile). 


rise was observed. During this temperature spike, the highest 
temperatures occurred at the trailing edge, particularly near 
the tip, where the airfoil was thinnest. 

At the second ignition spike evaluated, 1 .3 sec after startup 
(fig. 10(b)), the temperature gradient was spanwise, with the 
maximum temperature occurring at the leading edge near the 
airfoil base. This was expected because the high flow velocity 
yielded high heat transfer coefficients in this region. In 
addition, the large leading-edge offset caused a large gradient 
between the leading edge and the midspan surface. The 
temperature at the leading-edge base was approximately 
667 °C; during the first ignition spike, the temperature was 
approximately 870 °C. The trailing-edge region experienced 
much lower temperatures than it did during the first ignition 
spike. At 4.5 sec after startup (fig. 10(c)), a steady-state 
situation, or cruise, was encountered. Again, the hot-spot 
locations were at the leading and trailing edges at the airfoil 
base; the remaining surface experienced a relatively uniform 
temperature distribution. A gradient of about 100 deg C existed 
at the trailing edge between the base and the midspan. This 
gradient was due to cold hydrogen gas cooling the blade root 
region. Also, the surface temperatures at the edges reacted 
nearly identically and were closely coupled with the transient 
gas temperature. This was expected because the high gas-path 
pressures resulted in extremely high heat transfer coefficients 
at the leading and trailing edges. The advantages of controlling 
the gas temperature spike, or overshoot, were reflected in a 
reduced leading-edge temperature at the 1 .3-sec mark (fig. 9). 
The leading-edge temperature was reduced by a substantial 
100 to 600 deg C, and the difference between the leading edge 
and the midspan surface was reduced from 200 to 1 12 deg C. 

At cutoff conditions, the heat transfer results demonstrated 
that the maximum temperature was always at the leading-edge 
base (fig. 10(b)). The trailing edge was a secondary high- 
temperature region. The temperatures would have been much 
lower at the airfoil-shank junction if the temperature effects 
of the platform and the shank were not simulated by boundary 
conditions. Adding the platform would increase the heat 
dissipation by conduction through the platform surface. 

Another factor unaccounted for in the scaling procedure was 
the effect of boundary-layer history, or the blade heat 
capacitance in response to gas temperature. When gas 
temperature increases abruptly, the boundary layer cannot heat 
as rapidly as the main-stream gas; consequently, the boundary 
layer shields the blade surface from the full effect of increased 
gas temperature. To examine the magnitude of boundary-layer 
history effects, a boundary-layer analysis should be conducted 
for the peak of the first ignition spike. Lower temperatures 
than those presented would be expected. Because there was 
insufficient information to permit a boundary-layer analysis 
at that time, an effort is now under way to collect the 
information needed to conduct transient heat transfer finite- 
element analyses for a complete blade model, including the 
platform and shank regions. 
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TEMPERATURE, 



°c. 

(°F) 

1142 

(2087) 

1119 

(2046) 

1097 

(2006) 

1074 

(1965) 

1051 

(1924) 

1028 

(1883) 

1006 

(1843) 

933 

(1802) 

961 

(1761) 

938 

(1720) 

916 

(1680) 

893 

(1639) 

870 

(1598) 

847 

(1557) 

829 

(1517) 

802 

(1476) 

702 

(12%) 

687 

(1269) 

672 

(1241) 

656 

(1213) 

641 

(1186) 

626 

(1158) 

611 

(1131) 

595 

(1103) 

579 

(1075) 

564 

(1048) 

549 

(1020) 

534 

(993) 

518 

(%5) 

503 

(937) 

488 

(910) 

472 

(882) 

777 

(1430) 

765 

(1409) 

753 

(1388) 

742 

(1367) 

738 

(1346) 

718 

(1325) 

707 

(1304) 

695 

(1283) 

683 

(1262) 

672 

(1241) 

660 

(1220) 

648 

(1199) 

637 

(1178) 

625 

(1157) 

613 

(1136) 

602 

(1115) 


(a) At maximum spike (0.6 sec after startup), 
(b) At second spike (1.3 sec after startup), 
(c) During cruise. 

Figure 10.— Airfoil temperature distribution. 
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Summary of Results PH 

A three-dimensional, finite-element heat transfer analysis 
was conducted on the first-stage blades of the high-pressure 
fuel turbopump for the space shuttle main engine by using the 
MARC finite-element code. The modified STAN5 code was used 
to run a boundary-layer analysis that predicted the convective 
heat transfer coefficients. Predicted high-temperature locations 
on the blade airfoil were evaluated for a corrected gas temper- 
ature profile and compared with a reduced or actual 
thermocouple readout. All transient results presented were 
obtained by scaling steady-state heat transfer coefficients on 
the basis of transient flow and temperature. These results will 
be used as thermal input for the proposed structural analyses. 
The following general conclusions were drawn: 

1 . The thermal response predicted from the finite-element 
analysis showed that the leading- and trailing-edge bases of 
the airfoil are the hottest locations. This was expected because 
high heat transfer coefficients were present in these regions, 
and hardware evidence from mission flights as well as 
experiments has verified these results. 

2. The use of the uncorrected, rather than the corrected, 
gas temperature profile resulted in a 14 percent decrease in 
the airfoil metal temperature at the leading edge as well as 
a 50 percent decrease between the midspan surface and the 
leading-edge base. 

3. Temperature distributions showed a chordwise variation 
at the first ignition spike and a spanwise variation thereafter 
into the cycle. A uniform temperature distribution was 
dominant at most of the airfoil surface during cruise except 
near the base, where a mixture of cold and hot gas was present. 
The coldest spot was always at the airfoil base because the 
boundary conditions there were cooler. 
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